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ABSTRACT 

One of the effective mecfianisms of neutrino energy losses in red giants, pre- 
supernovae and in the cores of white dwarfs is the emission of neutrino- 
antineutrino pairs in the process of plasmon decay. In this paper, we nu- 
merically calculate the emissivity due to plasmon decay in a wide range of 
temperatures (lO'^ - 10") K and densities (2 x 10^ - 10^'*) g cm . Numerical 
results are approximated by convenient analytical expressions. We also cal- 
culate and approximate by analytical expressions the neutrino luminosity of 
white dwarfs due to plasmon decay, as a function of their mass and internal 
temperature. This neutrino luminosity depends on the chemical composition 
of white dwarfs only through the parameter fi^ (the net number of baryons 
per electron) and is the dominant neutrino luminosity in all white dwarfs at 
the neutrino cooling stage. 

Key words: stars: neutrino processes - red giants - presupernova - white 
dwarfs. 

1 INTRODUCTION 

It is well known that neutrino emission plays an important role in the evolution of red giants, 
presupernovae, white dwarfs, and neutron stars. Neutrinos appear in a number of reactions 
in dense stellar matter (see, e.g., Yakovlev et al. 2001) and freely escape from the star, 
producing a powerful mechanism of their cooling. One of the effective neutrino generation 
mechanisms is the plasmon decay. 

* E-mail: kantor@mail.ioffe.ru (EMK); gusakov@astro.ioffe.ru (MEG) 
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In contrast to ordinary photons in vacuum, plasmons, which are quanta of electromag- 
netic field in a plasma, can be not only transverse (in this case two polarization vectors of 
plasmon are perpendicular to wave vector), but longitudinal as well. The longitudinal plas- 
mons appear in the theory as a result of quantization of the well known Langmuir plasma 
waves. 

Plasmon can decay into a neutrino- antineutrino pair, 7 v-\-V. The appropriate neutrino 
emissivity was analyzed in a series of papers since 1963, when Adams, Rudcrman, and Woo 
had suggested this mechanism of energy losses in dense stellar matter for the first time. An 
account of these papers and references can be found in the review by Yakovlev et al. (2001) 
as well as in a recent paper by Odrzywolek (2007). Here we discuss in more detail only three 
papers which summarize and extend the results of previous works. 

Itoh et al. (1992) calculated the emissivity due to plasmon decay as a function of temper- 
ature and density and presented a table of numerical values and an approximate fit formula. 
Unfortunately, this approximate formula does not reproduce analytical asymptotes for the 
emissivity and thus can be applied only in a restricted region of temperatures and densities 
(near the maximum value of the emissivity). In addition, when calculating the emissivity, 
Itoh et al. (1992) used approximate expressions for the dielectric functions of electron gas 
and for plasmon dispersion relations which can be justified only at low enough temperatures 
(in a strongly degenerate electron gas). 

On the contrary, Braaten & Segel (1993) started with the most general expressions for the 
neutrino emissivity due to plasmon decay. They did not make any assumptions concerning 
degeneracy of the electron gas at calculating the dielectric functions and plasmon dispersion 
relations. To simphfy their analysis, Braaten & Segel (1993) suggested an elegant scheme 
to calculate approximately the dielectric functions, dispersion relations, and the neutrino 
emissivity. However, these authors did not present any tables with their numerical results 
or any approximate formula for the emissivity. Therefore, it is difficult to use their results 
in applications. 

Using the approximate method of Braaten & Segel, Haft et al. (1994) calculated the 
emissivity due to plasmon decay and fitted it by an analytical formula. This formula accu- 
rately describes the emissivity in a range of temperatures and densities where the plasmon 

decay is the most important neutrino emission mechanism. However, the fitting expression 
of Haft et al. (1994) does not satisfy the analytical asymptotes for the emissivity (they are 
presented in Section 2). 
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In this paper we would like to fill in the gaps in the literature devoted to the subject. 
We will {i) numerically calculate the neutrino emissivity due to plasmon decay making 
no assumptions concerning degeneracy or relativity of the electron gas; {ii) employ the 
approximate scheme of Braaten & Segel and find a fitting expression for the emissivity 
which reproduces the correct asymptotes. Thus, our main goal is to facilitate the use of the 
data on the neutrino emission due to plasmon decay. 

The paper is organized as follows. In Section 2 we present general equations describing 
the neutrino energy loss rate owing to plasmon decay. In Section 3 we give the fit expression 
for the plasma frequency which is a key parameter because the asymptotes of the emissivity 
depend on it. In Section 4 we present the fit expressions for the emissivity. In Section 5 
we apply the results of the preceding sections and find an analytical formula describing the 
neutrino luminosity of white dwarfs as a function of their mass and internal temperature. We 
summarize in Section 6. In Appendix A we present expressions for the dielectric functions 
of the electron-positron plasma. Finally, in Appendix B we describe a table of our numerical 
results. 

2 GENERAL EQUATIONS 

The neutrino emissivity due to plasmon decay can be presented as a sum of three compo- 
nents: the longitudinal component Q\ (due to decay of longitudinal plasmons); the transverse 
component Qt (the decay of transverse plasmons governed by the vector part of the weak 
interaction Hamiltonian) ; and the axial component Qa (the decay of transverse plasmons 
governed by the axial part of the weak interaction Hamiltonian). The component Qa is small 
and can be neglected (see, e.g., Kohyama et al. 1994). 

The emissivities Qt and Qi (per unit volume) are given in the form of integrals (see, e.g., 
Braaten 1991, Braaten & Segel 1993) 



Here, the integration is carried over the plasmon wave number k. In equations ([T]) and ([2]) 



Qo = [{m^cf/h^^] [G|/(967r^a)] (E,.^v) ~ 1-3858 x lO^^ ergs'^ cm"^; Gp = 1.436 x 10"' 



erg cm^ is the Fermi weak coupling constant; a = e^/{hc) ~ 1/137 is the fine structure 
constant; e and are the electron charge and mass, respectively; h is the Planck con- 




(1) 



(2) 
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stant; c is the speed of light; Cy ~ 0.9248 is the sum of squared normahzed vector 
constants Cy over all neutrino flavors. Furthermore, uJt{k) and uj\{k) are, respectively, the 
frequencies of transverse and longitudinal plasmons, which depend on the wave number fc; 
Zt(A;)~^ = d{uo1tt) / d{ijjl)] Zi{k)~^ = (ujf — k'^c^) de\/d{ujf), where et and e\ are the transverse 
and longitudinal dielectric functions of the electron-positron plasma, respectively. Finally, 
^B(t^t,i) = l/{^^P[^^t,i/(^B^)] — 1} is the Bose-Einstein distribution function for transverse 
or longitudinal plasmons; T is the temperature; ks is the Boltzmann constant; /^max is the 
maximum wave number at which the decay of longitudinal plasmon is still kinematically 
allowed by energy and momentum conservation laws. 

In the astrophysical literature the emissivity is presented as a function of temperature T 
and the effective mass density p, given by 

p = p/fie, (3) 

where p is the actual mass density; pe = J2i ^i^i/ ^i^i) is the net number of baryons 
per electron; Zi and Ai are, respectively, the charge and mass numbers of atomic nucleus 
species i; Ui is the number density of these species. Notice, that at densities higher than the 
neutron drip density pa ~ 4 x 10^^ g cm~^, free neutrons must be taken into account in the 
sum over i, in addition to atomic nuclei, when calculating p^. 
It is straightforward to verify that p can be rewritten as 

p ^ (ue - ne+) m^. (4) 

Here, rie and ne+ are the number densities of electrons and positrons; rrj-u is the atomic mass 
unit. 

The dependence of the emissivity Q = Qt + Q\ on p for temperatures T = 10'', 10^, 10^, 
and 10^° K is presented in Fig. 1. As seen from the figure, at fixed p the emissivity increases 
with the growth of T. If we fix T, the dependence Q{p) has a maximum. In the vicinity of 
the maximum the plasma frequency of the electron-positron plasma ujp is of the order of 
temperature, JiWp ~ k^T (see Section 5 for details). At high temperatures and low densities 
the emissivity ceases to depend on p (see equations El [IHl and [H] below) . In the figure this 
situation is illustrated by the upper curve, which is plotted for T = 10^° K. One sees that 
at p < 10*^ g cm~^ the curve tends to be horizontal. 

As follows from equations ([T]) and (I2l), for calculation of Qt and Q\ one needs to know 
the dispersion relations for transverse and longitudinal plasmons, uJi{k) and uJt{k), as well as 
the dielectric functions et{u!,k) and ei{u!,k). We calculated the dielectric functions et{uj,k) 



The neutrino emission due to plasmon decay and neutrino luminosity of white dwarfs 

30 




4 6 8 10 12 14 
Ig p (g cm-3) 

Figure 1. The emissivity Q = Qt + Qi versus p for T = lO'', 10*, lO'-*, and 10^° K. 

and e\{uj,k) for a wide range of densities and temperatures in the random phase approxi- 
mation and numerically obtained the dispersion relations and the plasma frequency uj^. The 
equations we used to compute the dielectric functions of the electron-positron plasma are 
given in Appendix A. These results were applied to calculate the integrals ([1]) and ([2]). In 
these calculations, we did not make any simplifying assumptions concerning the degree of 
degeneracy or relativity of the electron gas. The table with our numerical results can be 



found on the web, http : / /www. ioffe . ru /astro /NSG / plasmon /table . dat | (file table.dat). This 
table is described in Appendix B. 

The emissivities Q\ and Qt depend on two parameters characterizing stellar matter. For 
example, one may choose T and or T and p as proper parameters. Following previous 
results (see, e.g., Itoh et al. 1992), we take T and p as independent variables. It is convenient 
to introduce the notation / = hw^j {k-QT). 

The expression for the plasma frequency in the Braaten-Segel approximation has the 
form (see Braaten & Segel 1993) 



4a 



(5) 



where p, v = pc/E, and E = ^Jp^c? + m^c^ are, respectively, the momentum, dimensionless 
velocity, and energy of an electron or positron; up^E) = l/{exp[(_E' — p)/ {k^T)] + 1} is the 
Fermi-Dirac distribution for electrons; n-p{E) = l/{exp[{E + p)/{kBT)] + 1} is the Fermi- 
Dirac distribution for positrons; p is the electron chemical potential. 

In the region of relativistic temperatures {ksT ^ meC^) and under the condition k^T ^ 
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Pfc, the plasma frequency ([5]) has the asymptote 

2 _ Ana {k^Tf 

Here p-p = (STr^/i'^no)^/'^. For a degenerate electron gas pp is the usual Fermi momentum of 
the electrons. 

In the case when (i) the electron gas is degenerate (k^T <C ^/ppC^ + mld^ — rrieC^ and 
the contribution of positrons to Up can be neglected), or (m) the gas is non-degenerate, non- 
relativistic, and the temperature is not too high for the appearance of positrons [see, e.g., 
Landau & Lifshitz 1980, section 105], expression reduces to 

If the gas is non-relativistic (pp <^ iTiec), then this equation gives the well-known result, 
cjp = Ane'^ne/me. Notice that, since the contribution of the positrons to the asymptote (I7j) is 
negligible (ne+ <C Ue), Pf in this case can be approximately calculated as pf ~ [Sn'^h^ (rie — 
ne+)]^/^ = (Sn'^h^ p/rriuY^^ (see equation |4]). Introducing a new dimensionless parameter, 
Pf = {f^/rricc) (Stt^ p/rriu)^/'^, one can substitute {rrieCp-p) for p-p in the asymptote ([7]). 

Braaten & Segel (1993) developed a useful approximate method to calculate the emis- 
sivity due to plasmon decay. Below in this section we present some results obtained using 
this method (more details are given in the original paper of the authors). 

Using the method of Braaten & Segel, the emissivity can be expressed through the 
parameter f^,, which is a characteristic dimensionless velocity of electrons scaling from in 
the non-relativistic limit to 1 in the ultrarelativistic limit. 



Up 



(8) 

Here, the plasma frequency ojp is given by equation ([5]) while the frequency uj\ is 

^l = ^'^2 [ dp ^ l^^v' - v'^ [np{E) + ME)] . (9) 

In two limiting cases the neutrino emissivity due to decay of longitudinal and transverse 
plasmons can be calculated analytically. If the plasma frequency is much smaller than the 
temperature [/ = tujjp/{k-QT) <^ 1], then equations ([1]) and ([2]) can be simplified and writ- 
ten as 

Qt = go 4C3/3V^ (10) 

Qi = Qo A{v.)f. (11) 
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Here, — 1-202057 is a value of the Riemann zeta-function and the function (3{v^) equals 



3 

2^ 



1 



1 



In 



1/2 



(12) 



2v^ I - 

In the non-relativistic limit (f* 0) it reduces to /5 = 1, while in the ultra-relativistic limit 
{v^ — > 1) one has (3 = a/3/2. Furthermore, A{v^,) is a smooth function of f*, changing from 
8/105 ~ 0.076 at v^.. to 0.349 at f,,, ^ 1. If the plasma frequency is much greater than 
the temperature (/ ^ 1), then the integrals ([1]) and ([2]) can be taken analytically, 



Qt — Qo 

Qi = Qo 

where bi = 



knT 



rUpC^ 



hf'-^expi-f), 



hf'-'exp{-f), 



(13) 
(14) 



^{1 + ?;2/5)-3/2 and h = i3v^j5)-^/\ 



3 FIT FOR PLASMA FREQUENCY 

To simplify subsequent analysis we derived an analytical formula which approximates the 
plasma frequency ([S]) in a wide range of temperatures T = (10^ — 10^^) K and effective 
densities p = (2 x 10^ - 10^^) g cm . This range of parameters includes all possible limiting 
cases of degenerate, ultrarelativistic, as well as of non-degenerate non-relativistic electron 
gas. We calculated the emissivity on a dense grid of mesh points (with the steps 0.2 in IgT 
and Igp). The root mean-square relative error of our approximation is 0.4%. The maximum 
error of 1.4% is at IgT = 9.0 (K) and Igp = 2.4 (g cm~^). The fit reproduces the asymptotes 
from Section 2. The squared plasma frequency can be approximated as 

u^l = (^) ' V^asy^ + [asy, (1 - C D)]^. (15) 

Here, asy^^ = 4a/(37r) p| / a/ 1 + Pp is exactly the low-temperature asymptote (171) [we recall 
that pp = {h/mec) {3n'^p/m^^y/'^], while asy2 is given by 

-10 

(16) 

with t = k-QT / {meC^) . In the high-temperature limit, asy2 transforms into the asymptote 
([6]). The fit parameters pi, p2, and ps equal pi = 1.793, p2 = 0.0645, and ps = 0.433. 
The function C in equation (JT5l) is written as 



Ana 



asy2 



-P2 



t P2 
P2 



1 + 



P3 



it/^/p'2 



iPi 
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where 

1 +P7 (1 +J95 

C2 = P8+P9 (19) 

PlO + P 

with p4 = 0.01139, P5 = 2.484 x 10^ = -0.6195, p-j = 0.0009632, pg = 0.4372, pg = 1.614, 
and pio = 8.504 x 10^ 

At low temperatures the plasma frequency in the first approximation depends only on p 
and we have C = 1. The function D in equation f[T^ has the form 

^ = , (20) 

At high temperatures (when t ^ 1 and the electron gas is non-degenerate) the fit (fT5|) 
reproduces the high-temperature asymptote (jH]). At low temperatures (a degenerate gas or 
a non-degenerate non-relativistic gas; positrons can be neglected) the fit (IT3]) transforms 
into the analytical asymptote ([7]), which depends only on p. The function D is designed in 
such a way to reproduce not only the asymptote ([7]) of plasma frequency but also the first 
temperature corrections to Up. For a degenerate electron gas, the expansion parameter is 
k^T / p; for the non-degenerate non-relativistic gas it reduces to k-BT/irricC^). 



4 FIT FOR THE NEUTRINO EMISSIVITY 

In this section we present an analytical formula which approximates the results of numerical 
calculations of the emissivity Q = Qt + Q\ (per unit volume) and reproduces the asjTiiptotes 
from Section 2. The approximation was made in a range of temperatures T = (10^ — 10^^) K 
and effective densities p = (2 x 10^ — 10^^) g cm~^. The emissivity Q{p,T) was calculated 
on the same grid points as the plasma frequency (Section 3). At / = hw^/ik-oT) > 20 the 
accuracy of our fit is only logarithmic. However, in this case the emissivity Q is exponentially 
small, Q ~ exp(— /). 

The fit for the emissivity can be presented in the form 

Q = Qi + Qt = Qot' {Wt + Wi) exp(-/), (22) 
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where, as before, t = k^T / {nieC^) and we define 



Wt = asYti + asyt2 exp 



Q3 



(23) 



^ asyi2 [asyii + g4 (1 + gs vl^?'^ f^] 

asyi2 + [asyii + g4(l + g5t'PPr]' ^ ^ 



In equations fl23l) and ^ 
asyti = ai/^ asy^ = a2/^ (25) 
asyt2 = hf-\ asyi2 = 62/'-', (26) 



the functions 61 (t;*) and &2('i^*) are the same as in equations ( fT3il and ( 031) : the function Piy^) 
is given by equation (fT^ . At / ^ 1 equation (122]) transforms into 

Q = Qt = Qot'asyti = got' 4C3/3V' (28) 

(compare with the asymptotes [TUl and [TT]) . At / ^ 1 one has 

Q = got'(asy,2 + asy,2) exp(-/) = Q,t\h + &2)/'-' exp(-/) (29) 

(compare with the asymptotes [T3] and [T^ . 

When calculating the emissivity from equation fl22]) one should use the fit f[T^ for the 
plasma frequency cjp and the following fit for the characteristic velocity f*. 



1 + Si t*2p^3 



(30) 

where -up = Pf/a/1+Pp; Si = 9.079; S2 = 1.399; S3 = —0.06592. The root mean-square 
relative error of this approximate formula in the chosen range of T and p constitutes 1.4%. 
The maximum fit error is equal to 5.4% at IgT = 8.4 (K) and Igp = 3.8 (g cm~'^). 

In addition, it turns out to be necessary to use a special approximate formula for the 
function /5^(f*) from which the fitting expression ( l22l) depends on [a simple substitution of 
equation (|30l) into ( fT2l) and subsequent calculation of [3^ results in large errors], 



= f^'i^^) + [3.375 - (3%v^)] j^;-^j^y (31) 

where Vi = 0.3520; r2 = 1.195; = —0.1060. The root mean-square relative error of this fit 
constitutes 2.5%, the maximum error of 8.3% is at IgT = 9.6 (K) and Igp = 5.0 (g cm~'^). 
The function (I3T1) was approximated in the same temperature and density range as the 
parameter and the emissivity Q. 

The use of approximate formulae f|T5|) and fl30l) - fl3T|) leads to the following values of fitting 
parameters . . . ,(?5 (see equations [23] and [2l]) , minimizing root mean-square deviation of 
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Figure 2. Relative deviation <5 = (Qut — Qnum)/<9num versus p for T = 10^, 10^, 10^, and lO^'' K. Here, Qnum is the emissivity, 
numerically calculated in this paper. For Qut we take one of the emissivities obtained either from the fit formula ||22|I (solid 
lines); or from the table of Itoh et al. (1992) (long dashes); or from the fit formula of Itoh et al. (1992) (dots); or from the 
fit formula of Haft et al. (1994) (short dashes); or from the approximate formula from the review of Yakovlev et al. (2001) 
(dot-dashed lines). 



the emissivity, provided by equation fl22|) . from the numerical values, 
qi = 0.7886, = 0.2642, 

gs = 1.024, g4 = 0.07839, = 0.1784. (32) 

The root mean-square relative error of the approximate formula fl22l) with these coefficients 
is 4%, the maximum error is 7.9% at IgT = 8.4 and Igp = 6.4. 

In Fig. 2 we compare our numerical results for the emissivity Qnum with the results taken 
from the literature (corresponding emissivities are denoted as Qut)- The figure presents the 
relative deviation 6 = (Qht — Qnum)/Qnum as a function of p for a set of temperatures 
T = 10^, 10®, 10^, and 10^° K. The solid curves demonstrate relative deviations of the 
approximation (l22l) . suggested in this paper, from our numerical results Qnum; the long dashes 
show deviations from numerical calculations of Itoh et al. (1992) [taken from their table]; 
the dotted curves correspond to relative deviations calculated using an approximate formula, 
suggested by Itoh et al. (1992); the short dashes describe relative deviations calculated from 
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a fitting formula of Haft et al. (1994). Finally, by the dot-dashed curves we show relative 
deviations calculated from the approximate formula for the emissivity given in the review 
of Yakovlev et al. (2001). In that review it is recommended to use the formula only for 
p > 10^ g cm~^ and for strongly degenerate electrons. Prom the analysis of Fig. 2 a number 
of conclusions can be inferred: 

(i) The approximate formula obtained in this section is in good agreement with the 
results of numerical calculations as long as / < 20 (at greater /, that is at higher densities, 
the solid curve tends to go upward). 

(a) Our calculations agree with results of Itoh et al. (1992) in the range of parameters, 
where the electron gas is strongly degenerate and the emissivity is not small. However, as 
follows, for example, from Fig. 2 at T = 10^ K, some our results deviate from those of 
Itoh et al. (1992) for p ~ 10"*^^ g cm~^. For this case, matter is strongly degenerate so that 
the simplified assumptions, made by Itoh et al. at calculating the emissivity, could not lead 
to such deviations. (Let us note that Itoh et al. used the dielectric function, calculated by 
Jancovici 1962 for a strongly degenerate electron gas, see Appendix A.) Taking into account 
that our numerical results at such densities and T = 10^ K do not differ from the analytical 
asymptote for the emissivity by more than 10%, the results of Itoh et al. (1992) in the 
indicated parameter range seem less accurate than ours. 

(Hi) The fit formula of Itoh et al. (1992) satisfactorily describes the results of numerical 
calculations only near the maximum of the emissivity (when / ~ 1). 

(iv) The fit formula of Haft et al. (1994) agrees well with our numerical results in the 
same region of temperatures and densities in which the numerical results of Itoh et al. (1992) 
agrees with our numerical results. 

(v) The approximate formula from the review of Yakovlev et al. (2001) becomes inac- 
curate at high temperatures (i.e., T = 10^° K) and low densities (p ~ 10^ g cm~^). This 
approximate formula is valid only for strongly degenerate electrons, while the electron de- 
generacy becomes mild at high T and low p. 

Summarizing, as follows from Fig. 2, the results of various authors arc in satisfactory 
agreement in the ranges of T and p where the process of neutrino emission due to plasmon 
decay is the most efficient mechanism of energy losses in dense stellar matter. 
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5 THE NEUTRINO LUMINOSITY OF WHITE DWARFS 

Let us apply the results of Section 4 to analyze the neutrino luminosity of white dwarfs. As 
will be argued below, the neutrino luminosity due to plasmon decay only weakly depends 
on a specific model of a white dwarf. Thus, it can be considered as a universal function of 
the white dwarf mass M and its internal temperature T. Here we calculate this universal 
function and approximate it by a convenient analytical formula. 

As is well known, the thermal evolution of a white dwarf consists of two stages, the 
neutrino cooling stage (where cooling is mainly realized through the neutrino emission from 
the entire stellar body) and the photon stage (the main energy losses through the photon 
radiation from the stellar surface). A transition from one stage to the other occurs at the 
stellar age r ~ (10^ — 10*) yr, when the surface temperature of a star equals ~ 2.5 x 10^ K 
(for a hydrogen or helium atmosphere white dwarf, see, e.g., Winget et al. 2004). 

At the neutrino cooling stage the main mechanism of energy losses is the neutrino emis- 
sion due to plasmon decay (the second important process - the neutrino bremsstrahlung in 
collisions of electrons with atomic nuclei - is 10-100 times weaker, see Winget et al. 2004). 
We numerically calculated the neutrino luminosity Lj,{M,T) of white dwarfs caused by the 
decay of plasmons. When doing the calculation, wc made the following assumptions. First, 
to obtain the density profile inside a white dwarf we assumed that the pressure is fully 
determined by degenerate electrons. Second, the stellar core was assumed to be isothermal, 
which is a good approximation for not too young white dwarfs (r > 10 — 1000 yr). Third, 
we neglected beta-captures when calculating the structure and luminosity of massive white 
dwarfs. Beta-captures lead to softening of the equation of state, and influence the hydro- 
static structure of a star. In addition, they change stellar chemical composition, affect the 
number of nucleons per one electron, /le, and, consequently, the quantities p and L^,. How- 
ever, because the neutrino luminosity is the integral characteristic of a star, it should not 
strongly depend on these simplified assumptions. 

In Fig. 3 we present the neutrino luminosity L^, as a function of stellar core temperature 
T for white dwarfs with the masses M = OAMq, O.GMq, Mq, and 1.4M0. 

The results of numerical calculations of in the range of temperatures T = (3 x 10^ — 
5 X 10^) K and masses M — (0.4 — 1.3)M0 were approximated by the formula 



ke{l + A;5M)22/(3fc7)T|2/{3fer) ^ ^^^22/(3^7) 
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Figure 3. The neutrino luminosity Ly versus internal stellar temperature T for white dwarfs with M = 0.4Mq, 0.6Mq, Mq, 
and 1.4M0. 



where M = M/Mq, Tg = T/{10^K), and 
h = 1.050, k2 = 11.86, h = 5.901, 

h = 1.010, h = -0.5448, ke = 2.777, kj = 5.635. (34) 

For white dwarfs with M = (1.3 — 1.4)M0 the neutrino luminosity in the same range of 
temperatures T = (3 x 10'' — 5 x 10^) K is given by 

L.2(M,T) = 10^9- '-^^^ -J, ergs-\ (35) 

where 

h = 2.777, k = 25.13, k = 3.095, U = 7.585, h = 7.381. (36) 

The maximum error of the fit expressions (133!) and (135!) does not exceed 14%. Unfortunately, 
these two approximations do not match at M = 1.3M0. Thus, to calculate the neutrino 
emissivity for a white dwarf with the mass M G [I.28M0, 1.32M0], we recommend to use a 
linear interpolation 

L.3(M,T) = L.i(1.28M0,T) + ^^2(1.32M0, T)^- L.i(1.28M0, T) ~ _ ^ ^g), (37) 

This interpolation does not affect the maximum fit error which remains to be 14% at T = 
2.38 X 10* K and M = 1.34M0. 

As seen from equations (133|) and (!35|) . in the limit of high temperatures ~ T^, while in 
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the limit of low temperatures Ly ~ T^^^^. Let us demonstrate how to obtain this temperature 
dependence from simple physical arguments. 

At high temperatures, the internal stellar temperature T is much greater than the plasma 
frequency tUpo in the center of the star. Since the plasma frequency of degenerate electrons 
becomes smaller as the density decreases (see equation [7]), we have k-oT ^ hup throughout 
the star. In this case the neutrino emissivity of an arbitrary volume element in the star is 
given by asymptote ( l28l) and the luminosity equals 

L.^^^3 7^,QoT' [ P'u^ldV. (38) 

[TTIqC ) Jstar 

Here the integral is taken over the volume V of the star. Since the plasma frequency Up and 
the parameter /3 depend only on p (see equations [71 and [T2|) . one gets L^, oc T^. 

In the low-temperature limit, when k^T ^ fiUpQ, the main contribution to the lumi- 
nosity comes from a thin spherical layer of width h, in which hup ~ A^b^- This layer is 
situated in the outermost part of the stellar core, where the electrons form a degenerate, 
non-relativistic gas. Indeed, if we move from this layer to the stellar center, Up will increase 
while the emissivity will be exponentially suppressed, Q ~ exp^—hWp/kBT), in accordance 
with equation (1291) . If we move from the layer to the stellar surface then the emissivity will 
also decrease (see asymptote but in a power-law fashion, Q ~ P^ujp = (/3 = 1 for the 
non-relativistic electron gas, see equation [T2l) . Therefore, the emissivity will have a maxi- 
mum in a layer in which hup ~ /cb^, and the neutrino luminosity of a star can be estimated 
as 

U - 4^3/3^ Qo (^) ' (^) ' ^rcR'h ~ 16n^, Q, ' R^K (39) 

where R is the white dwarf radius. An order of magnitude estimate gives the characteristic 
width h of the layer, h ~ ci;p/(da;p/dr). Using the hydrostatic equilibrium equation and the 
scaling relations for the plasma frequency (see equation [7j) Up oc p^^^ and pressure P oc p^^^ of 
the degenerate non-relativistic gas, we get h oc p^/^ oc oop^^ oc 

T^/3_ Consequently, oc T^Vs^ 

in agreement with the estimate ([32]) • 

Let us note that the plasmon decay neutrino emissivity and hence the luminosity of the 
star depend on the effective density p, which is related to the real density p by equation Q, 
p/p = l/iUc = ^iZiiUi/ (YliiAini). In white dwarfs with any reasonable chemical composi- 
tion, the mass number Ai of atomic nuclei species i is always twice as much than the charge 
number Zi (recall that we neglect beta-captures). Thus, the ratio p/p is equal to 1/2. We 
used this ratio in all our calculations. 
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6 SUMMARY 

We have calculated the neutrino emissivity Q due to plasmon decay in an electron-positron 
plasma making no assumptions about degree of degeneracy or relativity of the electron gas. 

When calculating the emissivity one needs the plasma dielectric functions as well as the 
dispersion relations for transverse and longitudinal plasmons in a wide range of temperatures 
and densities. In particular, we have calculated the plasma frequency Up and fitted it by 
an analytical formula. This formula reproduces the main asymptotes for cUp (degenerate, 
ultrarelativistic or non-degenerate non-relativistic electrons, see Section 3). 

The results of numerical calculations of the neutrino emissivity were also approximated 
by a convenient analytical expression. It satisfies the asymptotes in various limiting cases 
(Section 4, also see the paper by Braaten & Segel 1993). The approximation is valid for 
T — (10^ — 10^^) K and p = (2 x 10^ — 10^^) g cm~^. The root mean-square relative error 
of the approximation does not exceed 4% for those temperatures and densities, for which 
/ = hwp/ (kBT) < 20 [while at / > 20 the emissivity is exponentially small, Q ~ exp(— /)]. 

The fit expression for the emissivity was used to calculate the neutrino luminosity of 
white dwarfs (Section 5). This neutrino luminosity was fitted by analytic formulas and 
presented as a function of white dwarf mass and its internal temperature. It is shown that 
the neutrino luminosity depends on the chemical composition of a white dwarf only through 
the parameter /i^ which is equal to 2 for reasonable white dwarf compositions. 

The results of this paper can be used in a number of applications, in particular, in 
modelling of the evolution of red giants or presupernovae as well as in the cooling theory of 
white dwarfs (see, e.g.. Haft et al. 1994, Winget et al. 2004). 
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APPENDIX A: DIELECTRIC FUNCTIONS OF ELECTRON-POSITRON 
PLASMA 

Using the density matrix formalism we calculated the dielectric function of the electron- 
positron gas in the first order of perturbation theory. The longitudinal e\{uj^ k) and transverse 
et(cu, k) components of the dielectric tensor can be written in the form [c = h= k^ = 1) 



Ana 



E 



d^p 1 
(27r)3 Ep+kEp Ep+k - Ep-io -iS 



■p+k 



2 ^ ,J +{p-k)+ Ep+kEp - E 



1 - 



A;2 
Ana 



(Al) 



E 



1 



d^jj 

(27r)3 Ep^kEp Ep+k - Ep-io -iS 



■p+k 



ipxkf 

k^ 



{p-k) + Ep+kEp - E^ 



(A2) 



Here, the summation is carried over electrons and positrons; Up = l/[exp(£^p =i= /i)/T + 1] is 
the Fermi-Dirac distribution function for electrons (in this case one have to choose the sign 
-) or positrons (the sign +); Ep = a/p^ + and -Ep+fc = a/ (p + ky + ml is the energy of 
an electron or a positron with the momentum p and p + k, respectively. 
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We have checked that equations (lAll) and (]A2p are equivalent to corresponding expres- 
sions for the dielectric function which can be obtained from the polarization tensor Tl^" of 
Braaten and Segel (1993) [see their equation Al]. 

The integration over the angles in equations (lAip and (lA2p can be done analytically. As 
a result, one obtains for real parts of e\ and ^t, 



^1 



a 



a 



dp [n^iEp) + ?2f(-Ep)] 
dpp^ [nY{Ep) +nY{Ep)] Rt, 



(A3) 
(A4) 



where 



Rt 



+ 



+ 
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Epk^ 



2Epk^p 

2 (u^ + e 
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Epk^ 



AEpk^p 



AEpk^p 



[k^ + Ak'^p^ + Ak'^EpUj - uj\2Ep + uf] In 
[k^ + Ak^p^ - Ak^EpUj - uj\2Ep - oof] In 



E'-{Ep + u) 



rp2 



(Ep 
(Ep 



El^u - {Ep 



UJ 



(A5) 



(A6) 



In equations flA3l) - (lA6P np{Ep) and np{Ep) are the Fermi-Dirac distribution functions for 
electrons and positrons, respectively; Ep±k = ^(p ± ky + ml is the energy of an electron 
or a positron with the absolute value of momentum equal to {p± k). 

Knowing the dielectric functions, the plasmon dispersion relations can be found from the 
equations 



ei{uj,k) = 0, u"^ et{uj,k) = k"^. 



(A7) 



If the electron gas is completely degenerate (T = 0), then the integrals in equations (]A3P 
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and flA4l) can be taken analytically. The result is 



a 

.2 , ,2 



+ 



+ 



1 (F - - 2ml) - + Ami 



3 - tu2) 

1 



ti;2 — fc2 



(A8) 
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2 1 3 p PFVp^ + m2 



-ie-uj^) sinh"!^ 
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A;3 
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24A;3 

The quantities L2 and L3 are 

-A;2 + _ 2A;pf)^ - ^oj'^ipl + '^e) 



In 
In 



;-A;2 + cu2)2 - A{uj^pl + ml + fcpp) 



(A9) 
(AlO) 

. (All) 

-A;2 + ti;2)2 — A{uJ^Jp1. + ml — kpp)^ 

The quantity Li depends on the sign of D = {00"^ — k'^){k'^ — 00'^ + 4ml). At D ^ one has 
-2mekpF + (A;2 + 2meUj — 00'^) ( \/pp + ml — 



Li = arctan 

+ arctan 

+ arctan 

+ arctan 
At D < 
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Note, that equations (]A8I) and (]A9|) for the dielectric functions agree with the well known 
results of Jancovici (1962) only at D < (see his equations Al and A4). At D ^ his 
expressions (Al) and (A4) are formally inapplicable (the real part of the dielectric functions 
in these equations becomes complex). In this case one should use our equations (lASP and 

dAlD. 

In addition, it may be useful to note, that the Jancovici's definition of the transverse 
dielectric function differs from a generally accepted one. His dielectric function e^^'i'^ is related 
to our dielectric function by sl^'^'^ = (fc^ — uj'^et)/{k'^ — uj"^)- 



APPENDIX B: DESCRIPTION OF A TABLE OF OUR NUMERICAL 
RESULTS 

The results of our numerical calculations are summarized in the table (file table.dat) which 



can be found on the web: http: / / www.ioffe.ru/astro/NSG/plasmon/table.datl 

The table consists of seven columns. In the first column, we present IgT (in kelvins); 
in the second column we give lg(p) = p//ic (g cm~^); in the third and fourth columns we 
present, respectively, the emissivities Qt and Q\ (erg s^^ cm^'^) due to decay of transverse 
and longitudinal plasmons; the fifth column is the plasma frequency tUp (s~^), which is 
numerically calculated from the exact dispersion relations ( ]A7[) [not using the Braaten- 
Segel approximation]; the sixth column is the same plasma frequency but calculated from 
equation ([5]) [the Braaten-Segel approximation]. Finally, in the seventh column we present 
the characteristic dimensionless velocity of electrons f^, = uji/ujp in units of c, calculated in 
the Braaten-Segel approximation (i.e. by making use of equations [5] and [9] for cUp and cui, 
respectively) . 
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